Center-of-mass corrections revisited: 
a many-body expansion approach 
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A many-body expansion for the computation of the charge form factor in the center-of- 
mass system is proposed. For convergence testing purposes, we apply our formahsm to the 
case of the harmonic oscillator shell model, where an exact solution exists. We also work out 
the details of the calculation involving realistic nuclear wave functions. Results obtained for 
the Argonne t;18 two-nucleon and Urbana-IX three-nucleon interactions are reported. No 
corrections due to the meson-exchange charge density are taken into account. 
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I. INTRODUCTION 

One of the successes of the shell- model picture has been 
the ability to calculate self-consistent densities for nuclear 
ground states that not only reproduce experimental bind- 
ing energies but also experimental charge radii of these 
nuclei and generally nuclear charge densities. The excel- 
lent agreement or remaining discrepancies have been a 
cornerstone for advancing our understanding of the nu- 
clear wave function. In particular, the ability to predict 
both, heavy and light nuclei, is taken as a confirmation 
of the quality of the effective nuclear interaction used in 
the calculations. For that reason it is useful to exam- 
ine the accuracy with which the nuclear densities can be 
calculated. 

For the proper description of the scattering process one 
assumes a nuclear wave function that factorizes into a nu- 
clear center-of-mass wave function, which is taken to be 
a plane wave, and an intrinsic wave function of coordi- 
nates relative to the center-of-mass. The difficulty lies in 
the ansatz of the wave function as a Slater determinant. 
Such a wave function generally does not factorize into 
a center-of mass wave function and a wave function for 
the nucleus relative to its center-of-mass. Furthermore, 
for the cases where it factorizes, the center-of-mass wave 
function is not a plane wave. While this is negligible 
for heavy nuclei, it is a significant correction for nuclei 
like 1^0. 

This problem has been known for a long time. It can be 
solved exactly for a single Slater determinant of harmonic 
oscillator single-particle wave functions. In that case it 
has been shown that the wave function factorizes with 
a center-of-mass wave function being a Gaussian. This 
allows us to calculate the form factor, i.e. the Fourier 
transform of the density, in the form 



Fsdiq) = 



Pint (q) 



(1.1) 



given in terms of the harmonic oscillator length parame- 
ter b. The calculation usually gives the form factor of the 
one-body density labeled Fgdiq) whereas the experiment 
requires the form factor with respect to the center-of- 
mass, labeled Fint{q)- Because of this exact result it has 
been customary to apply such a correction also in cases 
where the single particle wave functions are not harmonic 
oscillator wave functions and where the presence of cor- 
relations has been substituted by an effective interaction. 

An alternate way to deal with this is to calculate di- 
rectly the form factor in the center-of-mass system. This 
way the operator can be written as a series of one-body, 
two-body, to A-body terms. In this paper we first 
compare such an expansion with the exact result, for the 
case where such a result is available. We then apply the 
same expansion to a realistic wave function of ^^O [1] and 
compare it to the corrections implied by equation (1.1). 
This nuclear wave function was derived for ^^O using 
correlations of the form exp{S) together with the Ar- 
gonne 1)18 potential [2] that provides an excellent fit to 
the nucleon-nucleon scattering and thus must be consid- 
ered as a realistic interaction. Results corresponding to 
the inclusion of a phenomenological (Urbana-IX) three- 
nucleon interaction [3] are also reported. Thus, in this 
paper we hope to shed some light on the reliability of 
such center-of-mass corrections. 



II. THE FORM FACTOR OF THE DENSITY 

The charge form factor at momentum transfer q is 
jiven in Born approximation [4] by 



F..nt{q) = ('/'o I E 



(2.1) 
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where (pQ is the translationally invariant ground state, r' k 
the distance from the center-of-mass to the fcth "point" 
nucleon and fk{q^) the nucleon form factor, which takes 
into account the finite size of the nucleon k. 
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The ccntcr-of-mass correction has to do with the fact 
that the origin of the shell-model is not the same as 
the center-of-mass of the nucleus. Since the many-body 
Hamiltonian is not translationally invariant, then the 

model ground state ^g^^"* is not translationally invari- 
ant either, and thus can lead to incorrect description of 
observables, especially in small A nuclei. 

What we need to establish is the relationship between 
the model quantities expressed in terms of the coordi- 
nates of the laboratory system (r^., k ~ I ... A), and the 
intrinsic ones {r'k = Vk — Rem, k = 1 . . . A—1), measured 
from the center-of-mass of the nucleus 



Rr 



1 

7I: 



rfe . 



(2.2) 



fc=i 



Formally, this may be viewed as a change of coordinates, 
from the coordinates of the laboratory system fk to the 
coordinates of the center-of-mass system {Rcmj ^'fe}i fol- 
lowed by the removal of the dependence upon R^m from 
the model wave function ^q^^ , i.e. we have to construct 
the intrinsic wave function [5] 



i'i^\r'k) = J GiRem)^r\Rcm,r'k)dRem (2.3) 

independent of Rem, for an arbitrary function G{Rcm)- 
Note here that, in this formalism, the well-known 
Gartenhaus-Schwartz transformation [6,7] corresponds to 
taking G{Rcm) = 5{Rcm)- It is clear now that the arbi- 
trariness of the G{Rcm) function causes some troubles: 
Since there is no reason to choose a particular G{Rcm), 
it has been pointed out that the center-of-mass correc- 
tion for a given model wave function is not uniquely de- 
fined [5] . Nevertheless, the various recipes yield the same 
result in the limit of the exact wave function of a free nu- 
cleus [8]. 

The exact nuclear wave function $o consists of two 
factors, one of which is a plane wave in the center-of- 
mass coordinate, e^^'^""^, the other being the intrinsic 
wave function (/iq of the relative coordinates [9] r'k, 



{M), 



3'o(fi---fA) = e 



iPR„ 



Mr'i---r'A-i)- (2.4) 

For an approximate model wave function $q^^ however, 
all we can hope for is to be able to obtain the decompo- 
sition 

4"^^ = CpcmiRcm) <t>^^\r'i ■ ■ ■ r'A) , (2.5) 

which is approximately correct to the extent that the 
motion of the intrinsic coordinates and the center-of-mass 
are not correlated. Only then, the factorization 



Fsdiq) = Fcm{q) Fint{q) 



(2.6) 



is possible. To that approximation, and assuming that 
the model provides indeed a good description of the in- 
ternal structure of the nucleus (<l>o = ^o^^'' [10]), equa- 
tion (2.6) is valid with [8] 



k 

(2.7) 



and 



- /cb(^) I I 



(2.8) 



The form factor (2.7) can now be calculated directly 
by carrying out an expansion in terms of many-body op- 
erators 

k \ m^k I 

(2.9) 

Each exponential in equation (2.9) can be expressed in 
terms of the one-body operator which we define by 

f{q- Tm) = e-''-^^- - 1 . (2.10) 
With this we write the form factor as 
Fintiq) = Yl 

k 

X /e'«-'=(^-i)/^ H {l + nq-rm/A))\ (2.11) 

\ TOT^/C / 



or 



iq-fk(A-l)/A\ 



Fintiq) = E ^(«') - / 

k 

k m^k 

+ ^E^(«') E 

k m,n^k 

X (^e'f-^-'=(^-i)M/*(g . rmlA)r{q. Tn/A)) + ■ • ■ . (2.12) 

We intend to apply our formalism to the particular case 
of doubly magic nuclei (^®0). Thus, we can use the spher- 
ical symmetry of the nucleus to simplify calculations, in 
the sense that the form factor Fint{q) should be spheri- 
cally symmetric too, and we can in turn average the form 
factor over the directions of q. We introduce then 

F^^::\q)-^ jFWf)dn,. (2.13) 

This allows us to write the different terms in equa- 
tion (2.12) using the second quantization formalism, as 

follows 

1. one-body term 

E Uq') {a\0{q,n)\P) , (2.14) 

a/3 



with 



0{q,n) = Mqri{A-l)/A) 



(2.15) 
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2. two-body term 
^(2L + 1) 

L 

X X! •^«*^«^) |C'(5,ri,f2) I 75) aj^a^a^a^, 

(2.16) 



with 



0{q, ri, fa) = jL(<zri(A - I)/ A) fL{qr2/A) 

X (C^^^^ © C^"-^) . (2.17) 



3. three-body term 

J2 z^^-^^-^^{2L2 + 1)(2L3 + l)(i30L20 I LiO) 

X 1] /a(g') (a/37 1 0{q, n,f2, fa) I see) 
a/3759 c 



X aaa^a;fafa0a5 



with 



£'('7:^1,^2,^3) 

= JL, [qri {A - 1)1 A) fL, {qr2/A) U, {qr,/A) 



(ii) 



(2.19) 



and 



where we have introduced fi{qr) = ji{qr) — (5(o, 

ji{qr) and Cm = \J^^i^im{^) are the spherical Bessel 

functions of order / and the unnormahzed spherical 
harmonics of rank I and component m, respectively. 
Greek letters label the single-particle states |a) = 
\noc{laSoc)3ocmj^;Tam.rJ) , with s = i r 



and m-r = +\{—\) - for a proton (neutron). As a fi- 
nal remark, note that the conversion to second quanti- 
zation allows for all restrictions in the sums (2.12) to be 
dropped. 



III. HARMONIC OSCILLATOR SHELL-MODEL 
CALCULATION 

Equation (2.6) is always exact if ^o^^'' is expressed 
in terms of harmonic oscillator wave functions, provided 
that the centcr-of-mass wave function (f>cm is in one given 
harmonic oscillator state. Then, the extraction of the 
center-of-mass coordinate can be; done analytically. El- 
liott and Skyrme [11] have shown long time ago, that if 
the shell-model states are non-spurious, then the center- 
of-mass moves in its ground state and is described by the 
Is harmonic oscillator wave function 



{ A^ Y 
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(3.1) 



where h is the harmonic oscillator length parameter. The 
center-of-mass form factor can also be evaluated explic- 
itly 



Fcm{q) 



'/A 



(3.2) 



The correct translation-invariant form factor is thus given 
in terms of the shell- model form factor by 



Fint{q) 



FsM , 



(3.3) 



i.e. Fsd must be corrected by dividing through Fcm{q)- 
Note that, since the uniqueness of the procedure of car- 
rying out the center-of-mass corrections has been ques- 
tioned, the use of the equation (3.3) has been suggested 
even in the case of a more general nuclear structure 
model [9]. 

We exploit the analytical nature of these results by 
testing how fast does the many-body expansion (2.12) 
(2.18) converge. The shell- model wave function $0'*^'' for the 
harmonic oscillator potential is an independent particle 
wave function, represented by a simple Slater determi- 
nant of single-particle orbits. This state is what we 
shall call the uncorrelated ground state |0). By tak- 
ing the expectation value in the model ground state 



^0^^ = |0), of the one-, two- and three-body opera- 
tors in equations (2.14), (2.16) and (2.18), the following 
relevant expectation values are obtained 

{O\aiafs\O) = 6c,0 (3.4) 
(0 1 aj^a^a^a^, | 0) = SajSps - SasSpj (3.5) 
(0 1 ^t^b^isic^esis I 0) = 6aS {SpeS-yc - 5/3C^7f ) 

+ Sa<;{S0sSje - SffeS^s) ■ (3.6) 

Using these results and following a straight forward but 
laborious calculation, the translation-invariant form fac- 
tor for the harmonic oscillator shell-model can be com- 
puted completely up to the third-order in the many-body 
expansion (2.12). The various components involved are 
presented here, by their corresponding term of origin in 
the many-body expansion. Summations over all (nlj) 
indices are implicit. Notations are discussed in an Ap- 
pendix. 

a. One-body term. There is only one contribution to 
the one-body term of -F/^j"^ (q) 



(1)0 



ll01 = Uii{q')F:>: 



(3.7) 



where 7iOni{r) are the usual radial harmonic oscillator 
wave functions. Note that, in the previous equation, HOI 
is actually the Fourier transform of the one-body density 
folded with the appropriate nucleon form factor, i.e. 
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/•OO 

H01 = /p(g2) / p^P\r) joi^qr) dr 
Jo 

Jo 

(3.8) 

where p^^^ (r) and pp"' (r) arc the proton and neutron 
one-body densities, respectively, corresponding to the un- 
correlated ground state |0). 

h. Two-body term. Two components contribute to 
the two-body term of F^^^ '' (q) 

1. one component corresponding to the direct contrac- 
tion Sa-ySffS 

xn^^dr — Jnihjiyq )-^niliji,niliji-^n2l2j2,n2l2h ' 

(3.9) 

2. one component associated with the exchange con- 
traction SajS/^s 

H02ex = fnihji, 71212 j2iq^) X! ^'^^ + 

L 

f(2)L 



(1)L 

m/iji, "2(2^2 '^niliji,n2l2j2 



(3.10) 



where the pair of indices of the nucleon form fac- 
tor f{q^) indicate that the two orbits denoted as 
{nihji) and [7121232) have the same isospin. 

c. Three-body term. The three-body term contains 
six contributions to F^^^^ , out of which two are identical 
due to the fact that, in equation (2.18), the radial and 
angular parts of the operator dependent upon the coor- 
dinates of the 2nd nucleon are the same as the radial and 
angular parts of the operator dependent upon the coor- 
dinates of the 3rd nucleon. The different components of 
the three-body term (2.18) are listed below 

1. term 3.1 {SasS/Be^-ic) 



H03l — frnhjiiq"^) ^niliji,niiiji 

X 7(2)0 j{2)0 

n2l2j2,n2l2j2 na'ajajna'sja 



2. term 3.2 {SasSfi(;S-yg) 

H032 = -.fmhni<l'^)In^hjumhh 



(3.11) 



^(-1)^(2L +!)(/, 



(2)1, 



"2(2^2 ,"3(3^3 



f ; (3.12) 



3. term 3.3 (Sae^/ss^-yc) is equal to term 3.6 



HO33 = H036 

— —f ( r(2) 

— Jniiljl, "2(2^2 19 / Tl3i3j3,n3;3j3 

/ J ^ ' mil ji ,"2(2^2 mil ji ,"2(2^2 



; (3.13) 



4. term 3.4 (SaeS^^S^s) is equal to term, 3.5 

HO34 = HO35 = —fmhjt,n3l3j3{q^) 

X ^ (2L2 + 1) ^ (2L3 + 1) ^ i^+i^+La 



( L3 L 
\ 31 32 h ] 



j(l)L j(2)L3 J(2)L3 . 

nilijijUahja n2l2j2,niliji n3l3j3,n2l2h ' ^ '' 



In Fig. 1 we illustrate the convergence of the many- 
body expansion, for the case of the ^He and nu- 
clei, respectively. The solid line represents the exact 
form factor in the center-of-mass system, as given by 
Eq. (3.3). The agreement is excellent for a momentum 
transfer g < 3 fm~^, and remains reasonable good for 

up to 4 fm~^. It is expected that the size of the con- 
tributions due to correlations (as presented in the next 
section), is more important than the error made by ig- 
noring higher order terms in the many-body expansion 
(2.12). Also, it is worthwhile mentioning that a correc- 
tion expected to become increasingly important for high 
values of the momentum transfer, is the contribution due 
to the meson-exchange charge density [12]. However, the 
inclusion of this correction is beyond the purpose of the 
present discussion. 

We conclude that truncating the calculation at the 
third-order gives us a good approximation of the center- 
of-mass correction for the independent-particle model 
wave function case. Note that leaving out the three- 
body term in the case of the ^He nucleus, would result 
in an unacceptable description of the form factor distri- 
bution - false minima are located at a momentum trans- 
fer q as low as 3.6 fm~^ -, whereas in the case of the 
^^O nucleus, the charge form factor changes very little 
by including the three-body term. This is an indication 
that expression (3.3) can be viewed effectively, as a 1/A 
power expansion of the charge form factor. Therefore, as 
we consider the applicability of the expansion (3.3) for 
higher values of A, it appears that we can safely drop 
higher-order terms in the many-body expansion and still 
hope for a good description charge form factor. 

To conclude our study of the convergence of the many- 
body expansion, let us investigate the influence a given 
order of approximation has on the inferred mean square 
charge (rms) radius. It is well known that in the low q 
limit, the form factor may be be expanded in power series 
as 



Fint (q) 



1 
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q' ir') 



(3.15) 



and thus is a measure of the rms radius. Table (I) shows 
the convergence of the rms radius for the case of the ^He 
and ^^0 nuclei. These results show that the rms radius 
is little affected by any corrections beyond the two-body 
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term of the expansion (2.12). By including the three- 
body term in Eq. (2.12), the rms radius remains virtually 
the same in the ^He case, and changes by less than 1 % 
in the ^^O case. 



IV. REALISTIC NUCLEAR WAVE FUNCTION 
USING THE exp(5) METHOD 

We shall apply now our formalism to the case of a 
more complicated model wave function ^q^^ and the 
particular case of the ^^O nucleus. As advertised, the 
nuclear wave function ^Q*^-* = |0), has been obtained us- 
ing the coupled cluster method (or the exp(S) method) 
together with a realistic interaction [1]. The exact cor- 
related ground state ket wave function |0), is written in 
terms of the uncorrelated ground state |0), as 



|0)=e«^|0). 



(4.1) 



Here, is the cluster correlation operator, which may be 
decomposed in terms of p/i-creation operators (Oj = 1, 
0\ = at^a^i, Ol = at^at^a/^^a^i), as 



(4.2) 



n=0 



The expectation value of an arbitrary operator A in 
the energy eigenstate (4.1) may be written as 



A={0\e^Ae-^S^\0), 



(4.3) 



where similarly to S^, is defined by its decomposition 
in terms of p/ircreation operators 



(4.4) 



Therefore, the correct translation-invariant form fac- 
tor is given by the expectation value of the operator Fint 
in the correlated groimd state |0). As we have previ- 
ously [1] worked out the one- and two-body densities for 
the ground state, we can apply these results to evaluate 
the first two terms in this expansion. 

Using the definition of the one-body density 

p{r) = ^(0|5(f-f„)|0), (4.5) 

m 

together with the identity 

(4.6) 



dfe'i '^ 6{f-rk) , 
we can write the first term of Eq. (2.12) as 

Ai = /p(g2) J dfe'^-''^^-^^/^ p^P^r) 



Here, p^P\r) and p'^'^\r) are the proton and neutron 
ground state one-body densities, which include correc- 
tions due to 2p2/i, 3p3/i, and ApAh correlations. 

Similarly, we can write the second term as double in- 
tegral over the ground state two-body density, using 

p{n,r2) = ^ (0 I 5{n - rm) Kr2 - f„) | 0) . (4.8) 

mn 

Then, the second term of Eq. (2.12) becomes 
^2 = /p(g') j dr j dr' e'?-''^^-^^'^ Fiq-r'/A) 
p(P'P){f^f) + p(f.")(f,r') 
+ fn{q^) j dr j d? e'^-^^^-^y^ riq-r'/A) 



(4.9) 



(4.7) 



With these evaluations wc include all the terms that were 
included in evaluating the one- and two-body densities. 



V. RESULTS AND CONCLUSIONS 

The problem of center of mass corrections in calcu- 
lating observables has been worked out by expanding 
the center-of-mass correction as many-body operators. 
We have applied this expansion to the case of the har- 
monic oscillator where an exact solution exists. We found 
reasonable convergence in the case of harmonic oscilla- 
tor wave functions. Thus we have confidence that this 
method can be applied to general Hartree-Fock wave 
functions and in a situation where 2p2/i-correlations are 
present. 

Figures 2 and 3 show the various effects of the corre- 
lations on the internal charge form factor, corresponding 
to calculations using the Argonne wl8 with/without the 
Urbana-IX potential . Wc also compare the various ap- 
proximations of the form factor with the internal form 
factor suggested by Eq. (3.3), which in both cases is plot- 
ted as a dotted line. 

In the calculation of the translational invariant charge 
form factor correlations enter at two places. First, the 
calculation of the one-body operator (Al) includes ef- 
fects of all the correlations, because this term is sim- 
ply the Fourier transform of the one-body density. In 
Fig. 2, the solid and dashed lines represent the Fourier 
transform of the one-body density corresponding to the 
uncorrelated (|0)) and correlated (|0)) groimd state, re- 
spectively. These form factors are denoted SMl[po{r)] 
and SMl[p{r)]. Here, the main effect of the correlations 
is the shifting of the diffraction minimum by 5 % to the 
right. The new minimum is also predicted by Eq. (3.3), 
which also has a higher tail compared to SMl[po{r)] and 
SMl[p{r)]. 

Secondly, as any expectation value taken in the cor- 
related ground state, the center-of mass corrections are 
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modified due to the correlations. In Fig. 3, the solid 
and dashed lines represent the two-body approximations 
of the translational invariant form factor. Going be- 
yond the leading order {SM2) in evaluating the two-body 
term {A2), leaves the first diffraction minimum virtu- 
ally unchanged. However, the high q behaviour of the 
form factor, (g > 2.5fm~^), is dramatically affected. We 
can see that the Ai + A2 approximation of the inter- 
nal charge form factor exhibits a second diffraction min- 
imum, which has been observed experimentally by Sick 
and McCarthy [13] and its presence makes our theory 
credible. Physically speaking, the hole in the two-body 
density affects the center of mass motion and thus the 
center of mass correction to be applied. 



ACKNOWLEDGEMENTS 

This work was supported by the U.S. Department of 
Energy (DE-FG02-87ER-40371). Calculations were car- 
ried out on a HP-9000/735 workstation at the Research 
Computing Center, and a dual-processor 200 MHz Pen- 
tium Pro PC at the Nuclear Physics Group of the Uni- 
versity of New Hampshire. 



APPENDIX: NOTATIONS. 

We present here the various notations used in text. We 

have 

r(l)i 



= {'^hS{n^hjl), {7121202) + 1) 

/ HOnMr)nOn2i2{r)jL{qr{A - 1)/A)r^dr . 
Jo 



niliji,n2l2j2 
poo 

X 



'0 

,(2)1, 



nihil, v 



(2ji S, 



(miiii), ("2(2^2 



) + l) 



X / HOniii{r)nOn2i2{r)fL{qrlA)r''dr . 
Jo 



(Al) 



(A2) 



Here, the symbol S(niiiji),(n2i2j2) '-"^'^ when the set of 
indices (ni^iji) and (712^2^2) represent the same single- 
particle wave function, and zero otherwise. We also in- 
troduce 



f(1.2)L 
nihji,n2l2j2 

= {{hh)h\\c^'^\\{i2mCfiL2i 



2J2 



(A3) 



The reduced matrix element of the unnormalized spher- 
ical harmonic operator of rank k is 



(j{k) 



{hi)jb 



for 1^1 - ^2! < < ^1 + ^2 and \ji - 72] < k < ji +j2 , and 
zero otherwise. 
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q [fm '] q [fm '] 



(a) "He (a) Argonne vl8 



le+0 

le-1 ^ 

le-2 g- 

le-3 g- 

le-4 ^ 

le-5 ^ 
le-6 
le-7 

le-8 ^ 

le-9 ^ 

le-10 ^ 



2 3 

q [fm 



exact solution 
HOI 

H01+H02 
HOH-H02+H03 



1.5 2.0 2.5 

q [fm '] 



SMl[p„(r)] 
SMl[p(r)] 
Eq. (3.2.3) 




4.0 



(b) if^o 

FIG. 1. Convergence of the many-body expansion (2.12) of 
the charge form factor, for the harmonic oscillator shell model 



(b) Argonne d18 and Urbana-IX (preliminary) 

FIG. 2. ^''O nucleus: SMl[po{r)] and SMl[p{r)] form fac- 
tors compared with the internal form factor calculated ac- 
cording to Eq. (3.3). 
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TABLE I. Convergence of the mean square charge radius 
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q [fm '] 



(a) Argonne vl8 
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(b) Argonne d18 and Urbana IX (preliminary) 

FIG. 3. ^''O nucleus: Two-body approximations of the 
translational invariant form factor compared with the internal 
form factor calculated according to Eq. (3.3). 
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